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Recently, a systematic experiment measuring critical anomaly of viscosity of polymer solu- 
tions has been reported by H. Tanaka and his co-workers (Phys. Rev. E, 65, (2002), 021802). 

According to their experiments, the dynamic critical exponent of viscosity y c drastically de- 

, 

creases with increasing the molecular weight. In this article the kinetic coefficients renormalized 
by the non-linear hydrodynamic interaction are calculated by the mode coupling theory. We 
predict that the critical divergence of viscosity should be suppressed with increasing the molec- 
ular weight. The diffusion constant and the dynamic structure factor are also calculated. The 
present results explicitly show that the critical dynamics of polymer solutions should be affected 
by an extra spatio-temporal scale intrinsic to polymer solutions, and are consistent with the 
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experiment of Tanaka, et al. 
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§1. Introduction 

O 

Dynamic critical phenomena of classical fluids have been successfully investigated by the mode 
coupling theory 1-3 ) and the dynamic renormalization group theory. 4 ) Because of the divergence 
. of the correlation length £ of concentration fluctuations, the kinetic coefficients exhibit critical 
anomaly. According to the mode coupling theory and the dynamic renormalization group theory the 
anomalous kinetic coefficient for concentration behaves as e~ UXx , 3 ' 4 ^ where v is the critical exponent 
of £ nearly equal to 0.63 and e = (T — T c )/T c . Both theories predict that x\ = 0.95. Kawasaki 2 ) 
and Perl and Ferrell 5 ) predicted that the shear viscosity rj exhibits a logarithmic divergence: 



r\ = rj 



1 + 



(1) 



where r\ is the bare viscosity and Ao is a microscopic cut off wave number. Later, Ohta 6 ) pointed 
out that the viscosity anomaly is multiplicative and behaves as 

77 ~ e-"". (2) 
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The dynamic exponent y c has a universal number about 0.04. 7 ~ 10 ) One of the important notions 
in the studies of the static and dynamic critical phenomena is the universality. This insists that 
phenomena with the same set of critical exponents form a universality class. Dynamic critical 
phenomena of classical fluids belong to the model H universality in the Hohenberg-Halperin clas- 
sification. 4 ) 

It has been believed that the critical dynamics of complex fluids such as polymer solutions also 
belongs to the model H universality. Contrary to this conventional understanding, some experi- 
mental works have questioned that the dynamic critical phenomena of polymer solutions can be 
really categorized into the model H universality. n ~ 14 ) In fact the dynamic exponent y c obtained ex- 
perimentally for polymer solutions is smaller than that in classical fluids. n ~ 14 ) Recently H. Tanaka 
and his co-workers 15 ) have reported experimental work of dynamic critical phenomena of polymer 
solutions. Their accurate measurements for critical anomaly of viscosity show that the dynamic 
exponent y c significantly decreases as the molecular weight M w increases. That is, in polymer 
solutions with large molecular weight the critical divergence of viscosity is suppressed. They also 
measured the dynamic structure factor by the light scattering experiment. It has been found that 
in the case of large molecular weight the dynamic structure factor cannot be expressed by the 
Kawasaki scaling function even in the vicinity of the critical point. On the basis of such experimen- 
tal facts, they concluded that the critical dynamics of polymer solutions exhibits a non-universal 
nature and cannot be classified into the model H universality in a practical sense. 15 ) They sup- 
posed that such a nontrivial behavior is due to the dynamic coupling between critical concentration 
fluctuations and another slow viscoelastic mode which is intrinsic to polymer solutions. It should 
be noted that the viscoelastic effect and the break-down of the model H universality have been 
definitely confirmed in the studies for the phase separations of polymer solutions and dynamically 
asymmetric polymer blends. 16 ) 

Motivated by the experiments reported by Tanaka, et aZ, 15 ) we shall investigate the dynamic 
critical phenomena of polymer solutions by the mode coupling theory. In the present analysis we 
use the two-fluid model 17-19 ) as the basic equations describing the dynamics of polymer solutions. 
The two-fluid model was put forward by Brochard and de Gennes. 17,18 ) Later, the current form 
of the two-fluid model was derived by Doi and Onuki. 19 ) It has been successfully used to explain 
various phenomena such as shear induced concentration fluctuations, 19 ~ 22 ) non-exponential decay 
in dynamic scattering near equilibrium, 17 ~ 19 > 23 ) phase separation dynamics, 16 ' 24 ~ 26 ) and so on. 
However, the critical anomaly of polymer solutions has not been studied by the two-fluid model. 

Although the viscoelastic effect on the relaxation of the concentration fluctuations have been 
studied by a number of people within the framework of the two-fluid model (see for example 
ref. 19), the past studies were mainly performed in the systems where the mode coupling is not 
effective on the thermal relaxations, such as one phase region far from the critical point, very 
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viscous polymer solutions and polymer blends, and so on. Therefore, in these analyses the velocity 
fluctuations have not been explicitly taken into consideration. Consequently, they neglected the 
effect of the mode coupling on the thermal decay rates. 27 ) However, near the critical point in the 
semi-dilute polymer solutions, the dynamics of the critical fluctuations is strongly influenced by 
the non-linear hydrodynamic interaction arising from the streaming type mode coupling. As a 
result, the kinetic coefficients, observed near the critical point, are renormalized by the non-linear 
hydrodynamic interaction. In the present analysis such a non-linear hydrodynamic interaction is 
successfully taken into account by the mode coupling theory. We then show that the viscoelasticity 
affects the critical dynamics of polymer solutions, resulting in a suppression of the divergence of 
viscosity. This result properly explains the experiments reported by Tanaka, et a/. 15 ) 

The organization of the present paper is as follows: In §2 we briefly review the two-fluid model for 
the polymer solutions. In §3 we analyze the viscoelastic effect on the ralaxation of the concentration 
fluctuations and the velocity fluctuations within the linearized two-fluid model equations and derive 
the bare (unperturbed) propagators which give a starting point of the mode coupling analysis 
performed in §4. Here we shall evaluate the complex shear viscosity that is non-local both in space 
and time, which results from the dynamical asymmetry coupling between the viscoelastic stress 
and the velocity fluctuations. In §4 we investigate the critical anomaly, constructing self-consistent 
equations for the propagators by means of the mode coupling theory. The kinetic coefficients 
renormalized by the non-linear hydrodynamic interaction which is important in the vicinity of the 
critical point are calculated. In the final section we will give some remarks about the present 
analysis. 

§2. Basic equations 

In this article we use the two-fluid model for the basic equations to describe the dynamics of 
polymer solutions. Here we briefly survey the two-fluid model equations. The readers who want to 
know the details of the two-fluid model, see the original paper 19 ) or excellent reviews. 16, 26 ) 

Let v p (r, t) and v s (r, t) be the average velocities of polymer and solvent, respectively, and ip(r, t) is 
the volume fraction of polymer at point r and time t. We assume that the solution is incompressible 
and that polymer and solvent have the same specific volume. Thus ip satisfies a conservation law 
with v p : 

^V = -V-(Vv p ). (3) 
The volume-averaged velocity v = tpv p + (1 — ip)v s obeys the following hydrodynamic equation 

p^ = -Vp + r] V 2 v-ipV— + V-a p . (4) 

ot dip 

where p is the average mass density, r/o is the solvent viscosity, F is the free energy functional 
of the system, and a p is the viscoelastic stress. The hydrostatic pressure p is determined by the 
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incompressibility condition 

V • v = 0. (5) 

For slow motions, the two-fluid model gives the convective velocity of polymer v p as 

(I -ib) 2 SF 

where C(V') is the friction coefficient between the two components, and is given by 

C(V) = evrr/oC 2 - (7) 

Here is the blob length of order b/tp, with b being the monomer size. Eqs. (3) and (6) imply that 
the non-vanishing viscoelastic stress produces a diffusion. This effect is called the stress-diffusion 
coupling. The viscoelastic stress due to the network deformations is given by 

a'P = G(4>)W ■ (W-S). (8) 

This form is derived from the free energy functional introduced in eq.(10) shown below. Here G{ip) 
is the ^-dependent shear modulus, 5 denotes the unit matrix and W is interpreted as a long lived 
strain variable. We assume that the equation of motion of W is given by the following dynamic 
equation, 25,26 ) 

-QjT + (v p • V)W = (VvJ • W + W ■ Vv p ) - -{W - 5), (9) 
where r(ip) is the ^-dependent relaxation time of shear stress. 

§3. Linear Response 

In this section we examine the effect of the viscoelasticity on the relaxation of the concentra- 
tion fluctuations and the velocity fluctuations, based on the linearized two-fluid model equations. 
We then obtain the bare (unperturbed) propagators which give the starting point of the mode 
coupling analysis performed in the next section. Previously, Doi and Onuki 19 ) showed that the dif- 
fusion process is drastically influenced by the dynamical asymmetry coupling (a ^ 0) between the 
concentration fluctuations and the longitudinal mode of the viscoelastic force (stress-diffusion cou- 
pling). Their analysis predicted an additional wave number dependence of the diffusion coefficient 
due to the viscoelasticity, which was confirmed by the recent experimental studies (see for example 
ref. 27). However, in their analysis the very viscous systems without the velocity fluctuations were 
assumed, so the effect of the viscoelasticity on the relaxation of the velocity fluctuations (hydrody- 
namic relaxation) was not investigated. Here we take the transverse parts of the viscoelastic force 
and the velocity fluctuations into consideration. We shall then show that the dynamical asymmetry 
coupling between these transverse modes has a strong influence on the hydrodynamic relaxation, 
resulting in a wave number dependent shear viscosity. 
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The relaxation dynamics of the thermal fluctuations in a one phase region far from the critical 
point is still described by the linearized equations. Since the longitudinal modes do not couple to 
the transverse modes within the linearized equations, the diffusion process far from the critical point 
should not be influenced by the hydrodynamic relaxation. However, the situation is completely 
different in the vicinity of the critical point, since the non-linear hydrodynamic interaction, resulting 
from the streaming type mode coupling, becomes dominant. In the next section this will be shown 
by means of the mode coupling theory. 

Now we shall investigate the relaxation dynamics of concentration fluctuations and velocity fluc- 
tuations within the linearlized approximation of the set of equations represented in §2. We set 
To = r(V>o) and Go = G(ipo) with ipo being the average concentration. The free energy functional 
is assumed to be of the following Ginzburg-Landau type 20 ' 26 ) 

F{ V , W} = J dv[K oV 2 + ^col Wf + \Go5Wl], (10) 

where 5Wij = Wij — 5ij. The Fourier transform of an arbitary function g(r) is defined by 

<fr = J dve-^ r g{v). (11) 

In terms of the Fourier components, the linearized equations are given by 

p ~m = ~ Vok2 ^ + 1{5 " kk) ' k ' ^' (12) 

d 

— <5^ k = -r k 5Vk + aLk ■ k ■ <t£, (13) 
d 1 

+2G ar k fcjfcj<tyk - G a 2 Lk?(kifeio£ tlj + %^< H ), (14) 
<« = G 8W Kij , (15) 
where v k is the transverse part of v k , k = k/fc, and 

« = -f (16) 

Wo 

is the dynamical asymmetry parameter of the polymer solutions. 19 ' 24 ) The kinetic coefficient L is 
given by 

L = V 2 (l- Vo ) 2 /C( Vo ). (17) 
The decay rate in the absence of the viscoelastic coupling is given by 

r k = L(r + c £: 2 )fc 2 . (18) 

Here and in the following the Einstein summation convention has been adopted. The viscoelastic 
force f is defined as follows: 

f = V • a p . (19) 
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For the present purpose it is convenient to decompose the Fourier transform f k of the viscoelastic 
force into the longitudinal and the transverse components, 



f k = kk • f k + (6 - kk) • f k . 



(20) 



Let us set 



f» = kkf k , 
f£ =(J _ kk) . f k . 



(21) 
(22) 



From the linearized equations (12)-(15), we obtain the following two sets of equations of motion 
for the longitudinal modes and the transverse modes. 
For the longitudinal modes, we obtain 



3 7 
d~t Zk 



-T^ k -aLZ k + Rt(t), 



— (1 + 2^)Z k - 2G ar k ^^ k + K(t), 

TO 



where Z k = ik ■ and £„ e is a viscoelastic length or a magic length defined by 



£„ e = \J GoT a 2 L. 



(23) 
(24) 

(25) 



This length was first introduced by Brochard and de Gennes. 17,18 ) The physical significance of the 
viscoelastic length has been discussed by many researchers (see for example refs. 16 and 26). 
For the transverse modes, we obtain 



3yj 

dt 

dt k 



P P 
— (l + &tf)£-G tfv£ + R4(t). 

TO 



(26) 
(27) 



Here R z , R£ and R[ are random forces. 
First we formally solve eqs. (23) and (24) as 

SMt) = of (t)# k (0) + Of (t)Z k (0) 

+ f ' dsG^(t-s)Rt(s)+ f dsGl z {t-s)Rl{s), 
Jo Jo 

ZUt) = G^(t)# k (0) + Gf(t)2 k (0) 

+ f t dsG^(t-8)Ht(s)+ f dsG zz (t- s)R z {s), 
Jo Jo 

where the propagators G(i)'s are written as 



(28) 
(29) 



LL> + — UJ^ 



u+ - -(1 + 2f ve k 2 )} exp(-LU + t) + \-u>- + -(1 + 2^ e k 2 )} exp(-o;_t) 1, 

TO 1 TQ J I 
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Gf* z {t) 
Gl z (t) = 



aL 



UJ + — Ld- 

2G aT k k 2 

UJ + — LO- 
1 



exp(— oj+t) — exp(— uj-t) 

exp(— uJ+t) — exp(— uj-t) 
(uj + - T k ) exp(-u; + t) + {—uj- + r k ) exp(-w_i) 



iO + — UJ- 

The relaxation rates, uj + and are given by 



(30) 
(31) 

(32) 

(33) 



r k + l(i + 2^ 2 ) 

to 



± 



TO J 



4^ 

TO 



(34) 



The propagator G^(t) is the dynamic structure factor, and it is expressed by a superposition of two 
exponential functions. This form of the dynamic structure factor was first proposed by Brochard 
and de Gennes 17 ' 18 ) on the basis of a phenomenological argument. Afterwards, Doi and Onuki 19 ) 
derived the same form of the structure factor by analysing their two-fluid model equations which 
is employed in the present analysis (see also ref. 24). 
Similarly as above, we can formally solve eqs. (26) and (27) and obtain 



+ / dsHZ v (t - s)RJ(s) + / dsH^ (t - s)Rl(s), 
Jo Jo 

+ f dsH^{t - s)Rl(s) + / dsH^(t - s)R{(s). 
Jo Jo 

where the propagators il(t)'s are given by 



(35) 
(36) 



n + - 1(1 + £ e k 2 )] exp(-fi+t) + [-«- + -(1 + Ck 2 )] exp(-n_t)) 1, 



TO 



TO 



G n k 2 



i 



exp(— 17 + i) — exp(— 0_i) 



exp(— Q + t) — exp(— 



(fi _ 1?°*!) exp(-^ + t) + + ^) exp(-ft„t) 

P P 



n + - n 

The relaxation rates fi + and f2_ are given by 



(37) 
(38) 

(39) 

(40) 



P T 



r] k 2 



1 



P T 



AG k 2 



(41) 
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The viscoelasticity is characterized by the complex shear modulus G^_(lo) = G' k (uj) + iG^(uj), 
where G' k (uj) and G k (u;) are called the storage modulus and the loss modulus, respectively. 31 ) 
From the Eqs.(26) and (27), the total shear stress a^oo) is expressed as 

oiMiM = [ikiv£j(u) + ikjv^u)] , (42) 

where we have introduced the Laplace transforms of v k (t) and &k(t) as 

v^H = / dtwi{t)e~^\ (43) 
J o 

POO 

a k (cj) = / dta k (t)e- iu,t . (44) 
J o 

In the present system G^(cj) and G k (u) are given by 

G ' kH = ^ + e ve k*r (45) 

G£( W ) = u, [ m + Goro^^^^] , (46) 
which are non-local both in space and time. The complex shear viscosity is given by 

tfrM = (47) 

ILU 

The imaginary part of G k (u>) gives the shear viscosity for long time scale motion (u> = 0) as 



where we have defined r\ v = GqTq. The asymptotic behavior of %(()) is represented by 

(50) 



V*(Q) = Vo + Vp (^e«l), (49) 
~ + rip72TT2 ( fe ^e > !)• 



For polymer solutions the two-fluid model gives the viscoelastic length as 19 ' 26 ) 

U = {VpMkb- (51) 

Thus, Eq.(50) is also expressed as 

%(0)=%(l + ^), (*&e»l), (52) 

which explicitly shows that the hydrodynamic interaction becomes weak beyond the blob length 
It is worth mentioning that this property is due to the dynamical asymmetry coupling (a 7^ 0) 
between the viscoelastic stress and the velocity fluctuations. For polymer solutions, it is well known 
that the diffusion process is drastically influenced by the dynamical asymmetry coupling between 
the longitudinal modes, namely the so-called stress-diffusion coupling. 17 ~ 19 ) The present analysis 
shows that the dynamical asymmetry coupling between the transverse modes is also important for 
the hydrodynamic relaxation process. 
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§4. Mode coupling approach 

In this section we discuss the critical anomaly of the kinetic coefficients. If we take the whole 
non-linearity into consideration, many terms appear even in the second order in the perturbation 
expansion in terms of fluctuations and the analytical approach becomes very complicated. However, 
it turns out that there are only a few non-linear terms which crucially affect the critical dynamics. 
From the equations of motion represented in §2, we find that the non-linear terms, associated 
with the viscoelasticity, are so much smaller than the linear terms which are represented in the 
linearlized equations (12)-(15) that we may safely discard the non-linear terms arising from the 
viscoelasticity, and retain only two streaming type mode coupling terms which arise from the 
hydrodynamic interaction: One is the convection term of v in the equation for ip and the other is 
the osmotic pressure gradient term in the hydrodynamic equation. We consider these nonlinearity 
up to the second order of the perturbation expansion. As a result the equations of motion that we 
must analyze are given as follows: 

dv 1 - r f i 

= -r] k 2 v£ + i(5 - kk) ■ k • <r£ - / p(r + c p 2 )Y> k -p^p , (53) 

L Jp J 



r at 
a_ 

dt 



TO 

+2G aT k k i kj8^] i - G a 2 'Lk 2 {hh^ij + kjk[a^ H ), (55) 
<« = G 6W Kij . (56) 

As shown in the following sub-sections, our analysis based on the above set of equations well 
explains the experimental results. 15 ) 

4-1 Self- consistent Equations for Propagators 

In this sub-section we construct the self-consistent equations for the propagators within the mode 
coupling theory. 1 ' 2 ) By the simplest approximation we obtain the following set of self-consistent 
equations for the propagators: 

j^G k (t) = -i£ ■ G k (t) 

W r» /\ 2ri fA *2,<ltfk-p |2 



p 



J q ds jf k 2 [l - (p • k) 2 ]^^H; v (s)Gtt p (s)E ■ G k (t - s), (57) 



^H k (t) = -l£ ■ H k (t) 



f ds j p 2 [l-(p-k) 2 ] x 
Jo Jp 



k B T 

(IV'k-p| 2 )(IV'p| 2 )(7r^y - (| ^ |2) ) 2 cy( a )Gg p ( 8 )^ • HUt - s), (58) 
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where (• • •) denotes the equilibrium average, 



-G 



L 



■4>z 



I** 



Tk ah 



fH _ 
^k — 



and the matrix E is given by 







-I 


[ i 


% ) 





G k 2 ±(i+e ve k 2 ) 



E 



1 




(59) 
(60) 

(61) 



Equations (57) and (58) correspond to the diagrammatic one-loop approximation, neglecting vertex 
corrections. 

4-2 Renormalized Diffusion Constant 

First we calculate the diffusion constant and the dynamic structure factor in the vicinity of the 
critical point. The equations for the propagators can be written as 



iu5 + L£ + Vg(w)E] ■ G k (uj) = 5, 
iu5 + L£ + (lo)e\ ■ H k (uj) = 6, 



where 



G k (w) 



dte~ iujt G k (t), 



dte- iujt H k (t). 



Here the self-energies Ti^(uj) and Sj^(w) are given by 



4p Jo j p 

2\/l„;. |2 



(|Vk-p| 2 )(|^p|^)( 



(IV'kl 2 ) 

p 2 [l-(p-k) 2 ] X 

1 1 



2G^ (a)G e (a)e - 



(62) 
(63) 

(64) 
(65) 

(66) 
(67) 



(IV'pI 2 ) (l^k- P | : 

These equations are so complicated that it is hopeless to solve them analytically in general situ- 
ations. However, near the critical point we are allowed to investigate eqs. (66) and (67) by the 
following iterative manner. 

First we evaluate (w) assuming that the shear viscosity does not exhibit any critical anomaly. 
As we shall see in the following analysis, the decay of Gf^{f) is governed by the slowest mode near 
the critical point. Therefore, Hp V (s) in (66) is given by (37), so that we may set G^_ (s) = 1. 
Considering the decay of concentration fluctuation, to in (66) is much smaller than fi + or f2_, 
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whence we use the Markov approximation to set u> = in (66). That is, our present analysis can be 
regarded as the zeroth order approximation for the additional frequency dependencies of the decay 
rates. If we use the Ornstein-Zernike form for (IV'kl 2 ) oc (k 2 + £~ 2 ) -1 , the resultant equations can 
be integrated appropriately as follows: 



e£(o) 



l-(k-p) 



p V 

k B Tk 2 

p (Vo + V P )P 2 
k B Tk 2 

(vo + v p )p 2 



l-(k.p) 2 

l-(k-p) 2 



V0 + Vp + W&eP 2 1 + £ 2 ( k - P) 2 

i + ek 2 



i 1 + £ 2 (k - p) 2 

Vp^veP 2 



2 1,2 



l + ^k 



Vo + V P + VotleP 2 1 + ? 2 ( k - P) 



k B T 



K(x) + ^F(x,A) 
Vo 



(68) 



67r(?7o + v P )C 3 

where x = £k, k B is the Boltzmann constant, and the parameter A is given by 

A = -Ljl + ^. (69) 

£ve V VO 

By using eq.(51), in the case rj p 3> Vo A is practically the ratio of the correlation length £ to the 
blob length 



A 
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In eq.(68) the first term is the Kawasaki scaling function, given by 

3 



K(x) 



1 + x 2 + (x 3 - x _1 )tan 1 x 



(70) 



(71) 



In the present case an additional term F(x, A) arises from the viscoelasticity, which is given by 

a 2 -l (1 + x 2 ) 2 -^ 4 



+ 



2x 3 



AA 2 x 3 



tan x 



1 + x 2 - A 2 
2Ax 2 



2\2 



1 + 



:i+x 2 ~A 2 ) 

AA 2 x 2 



2Ax 



1 + x 2 - A 2 



tan 1 



2Ax 



1 + x 2 - A 2 



(1 + x 2 ) 2 - A A 
AA 2 x(l + x 2 ) 



1+ (l±x 



A 



2\2 



AA 2 x 2 



:tan _ lv ^TT) x 



where 



s = 



(1 + x 2 - A 



2n2 



4A 2 (l + x 2 ) ' 

Though the above expression is very complicated, its asymptotic form becomes simpler: 



££(0) - 



(72) 



(73) 



(74) 
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16£ 3 (% + r/ p ) 



A r 1 X 

3vr?7o A I £ £ b 



-x 



(k » -r). 



(75) 
(76) 



16£ 3 7? ~ & 

Here we have neglected the terms irrelevant to the following analysis. We plot the reduced self 
energy and its asymptotic form for various values of r] p /r]Q and A in Figs. l(a)-(f). We note that 
the dynamic scaling is violated and the self energy cannot be expressed by the Kawasaki scaling 
function. This is evident from Fig. 1(a) and Fig. 1(d) where the viscoelasticity is found to be 
relatively strong. 
The dynamic structure factor is given by 



u + — w- 



+ 2f ve k 2 )] exp(-cD + t) + f-w_ + -(1 + 2f ve k 2 ) exp(-u;_t) 



TO 



TO 



(77) 



where the renormalized relaxation rates, uj + and lD_ are obtained as 



± 



r k + s^(o) + -(i + 2^ 

TO 



\ 



r k + sG(o) + -(i + 2efe 2 ) 

TO 



4l 

TO 



r k + E£(o)(i + 2&Ar») 



(78) 



A comment will be made on the above results. Although the functional form of the dynamic 
structure factor G^(t) is the same as eq.(30) derived in the previous section, the relaxation rates 
Cj± of the critical concentration fluctuations obtained here are the ones renormalized by the non- 
linear hydrodynamic interaction. 
In the vicinity of the critical point £ k (0) dominates r k , so the characteristic frequencies become 



1 



u+ = -{l + 2e ve k z ), 

TO 

u_ - E«(0). 



The resultant structure factor is approximated as 

Gf (t)<*exp[-E£(0)t]. 



(79) 
(80) 

(81) 



That is, the critical mode associated with the concentration fluctuations can be free from the 
coupling with the viscoelasticity very near the critical point. Then, using eqs.(51), (70) and (74), 
the long-wavelength expression of the diffusion constant is given by 

k B T 



D 



1 + 



GntiVp + Vo) V VoA 



12 



^{v P + Vo) 

k B T 
^{V P + Vo) 



,V0 I 



6 



1+ % 2 % 



Vo 



(82) 



We point out that the temperature regime where the hydrodynamic interaction dominates the 
diffusion process (£^(0) ^ ^k) becomes narrower as the molecular weight increases, since r\ v oc 
M 3 31 ) 

4-3 Renormalized Viscosity 

In the previous sub-section, we assumed that the shear viscosity does not show the critical 
anomaly. However, as is well known, the shear viscosity exhibits a weak divergence at the critical 
point. In this sub-section, we investigate the critical anomaly of the shear viscosity. We can evaluate 
the self energy £^(0) by putting oo = in (67): 

P 2 [i-(p-k) 2 ]x 
i i 



4p Jo j p 

IV>k-pl 2 )(l^ P | 2 >(- 



-) 2 G^(s)Git p (s). 



(83) 



'(I^pI 2 ) (l^k-pl 2 )' 

Retaining the lowest power of k, the resultant integration can be performed analytically, using the 
results obtained in the previous sub-section: 

k B T£ 4 f A o . rP 1 



rS£(0) 



/■A pO 
fc 2"JcV"y- 307r 2 J rfP (T^2 p 2)2 E G( )> 

where Ao is the microscopic cut-off wave number. The above integrand can be estimated as 

W p 6 1 ^ VO + Vp ,5, , 1, 

30^ (1 +^V) 2 SG(0) 5vr(l + ^ A y P ' [P ^ 

^ 8(t/q + rip) 1 



8r? 1 , NN 1, 



where r = 8r/ p ^/37r?7o^4- In the case r\ v » r/o eqs.(51) and (70) give 

^/*Vt ~ ( Vv^ 



6. 



The anomalous part of the shear viscosity Ar/ is given by 

Thus, we obtain the total shear viscosity in the following form, 
rj = fj + Ar/ 
= fj\l + 



(84) 

(85) 

(86) 
(87) 

(88) 
(89) 



15vr 2 



3vr 



1 



401 + ^4 



2r + ln(|±^) + 5ln(6Ao) 



r/ 



(90) 
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where fj = r/ + rj p is the bare shear viscosity. In Fig. 2 we plot the reduced shear viscosity At] /fj = 
T]/fj — 1 as a function of rf p /r)o- It is evident that the critical divergence of the shear viscosity is 
suppressed as the molecular weight increases. 

For T] v 3> 770 and £„ e = £, it is readily shown that the resultant renormalized shear viscosity does 
not exhibit critical anomaly. 



77 

fj 



1. 



(91) 



For ?7 P = 770 and £ ^ 6ue> where the viscoelasticity is relatively weak, we obtain the anomalous 
shear viscosity as follows: 



77 

fj 



^ 1 + 



15vH 



In 



£A C 



(92) 



Because of the small coefficient of the logarithmic term, it may well be exponentiated with the 
small exponent x v = 8/157T 2 = 0.054 as 



V 
V 



£A 



VP 

(^Ao)T 



(93) 



In the low molecular weight limit (rj p = 0) the expressions reduce to those of the simple classical 
fluids, with y c = x v i> ^ 0.034. 32 ) 

The above results are due to the existence of an extra length-scale £„ e intrinsic to the viscoelas- 
ticity. As is well known, 1 ' 4 ) in the case of classical fluids the behavior of the critical anomaly of the 
shear viscosity does not depend on the detail of the material. That is, the behavior of the critical 
anomaly does not depend on the material parameter, and is universal (model H). Contrary to this, 
as shown in the present analysis, in the case of polymer solutions the critical divergence of the shear 
viscosity strongly depend on the molecular weight M w . Now, a few comments must be made on 
those points. When the correlation length £ is very large compared with the other length scale, it is 
shown that our results certainly exhibit the universal behavior. However, as noted by Tanaka, 15 ) it 
is difficult to experimentally access such a temperature region for polymer solutions whose polymer 
has a very high molecular weight. In this sense we can say that the model H universality does not 
hold, at least in a practical sense. ^ 

§5. Conclusions and Remarks 

In this article we have investigated the dynamic critical phenomena of polymer solutions. Based 
on the mode coupling theory, we derived a set of self-consistent equations for the self energies, taking 
the non-linear hydrodynamic interaction into account. Renormalizing the non-linear hydrodynamic 
interaction terms which are important in the vicinity of the critical point, the transport coefficients 
were calculated. 



14 



Our analysis predicts that the critical divergence of the shear viscosity is suppressed as the 
molecular weight increases. This effect was indeed observed experimentally by Tanaka, et al. 15 ^ 
At the present stage the quantitative comparison with the experiments is difficult in the following 
reasons. There are many unknown parameters (for instance, rj p and etc.) which cannot be 
evaluated explicitly. In addition, in their paper 15 ) they introduced an effective exponent y c and 
then discussed the singularity of the shear viscosity. However, this does not relate directly to our 
result, eq. (90). 

We have also calculated the diffusion constant eq.(82) and the dynamic structure factor G^(t) 
with the renormalized relaxation rates u>±, eq.(78), by the non-linear hydrodynamic interaction. 
Those are new results for critical polymer solutions. I expect that the present theoretical analysis 
would stimulate further experimental studies. 

In the present analysis we have used the Markov approximation in order to evaluate the self- 
energies. Previously several authors investigated the memory effects in the critical dynamics of 
classical fluids and predicted that the effect has small contributions to the kinetic coefficients. 33 ' 34 ) 
These theoretical predictions have been successfully confirmed by the experiments. 35 ) In the case 
of polymer solutions the separation of the time scale between concentration fluctuations and the 
velocity is not clear as in the case of classical fluids. Hence there is a possibility that the viscoelas- 
ticity amplifies the frequency dependence of the kinetic coefficients. We will present more detailed 
studies for this elsewhere. 36 ) 

Finally, we note the following. It is well known that the dynamical asymmetry coupling (a ^ 0) 
between the concentration fluctuations and the viscoelastic stress gives the so-called stress-diffusion 
coupling. 17 ' 19 ^ In addition to this effect, by the present linear response analysis (§3), it has become 
apparent that the hydrodynamic relaxation process is also drastically influenced by the dynamical 
asymmetry coupling between the velocity field and the viscoelastic stress. The calculated bare 
shear viscosity has a wave number dependence as seen from eq.(48). 
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Figure Captions 

Fig.l. Reduced self energy Q(x, A, r/ p /ryo) = [K(x) + j^F(x, A)]/x 2 fromeq.(68) (solid line). The 
asymptotic form (broken line) is also plotted. We note that the dynamic scaling is broken down 
and the self energy cannot be expressed by the Kawasaki scaling function (dotted line), that is 
definite in (a) and (d) where the viscoelasticity is relatively strong. 

Fig.2. Reduced shear viscosity Arj/fj = rj/fj — 1 vs %,/r/o for various values of A. Here we set 
Ao6> = 10 and use eq.(70). It is evident that the divergence of the shear viscosity is suppressed as 
the molecular weight is increased. 
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